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determining the nuclear matrix elements for neutrinoless double beta decay. 
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1. Introduction 

With the discovery of neutrino oscillations, it has become evident that there are mass 
differences between the different eigenstates of neutrinos. This implies that neutrinos 
have a finite (though currently unknown) rest mass [U [2], where the upper limit is 
currently 0.28 eV/c^ f3j. Majorana had suggested 75 years ago [1] that the statistics for 
neutral, spin-1/2 particles might not be the same as those for charged ones, and allowing 
the possibility that neutrinos and antineutrinos could be one and the same. As long 
as neutrinos had zero mass, chirality distinguished neutrinos from antineutrinos. But 
with a finite rest mass established, the handedness of neutrinos is no longer an intrinsic 
property, making Majorana's hypothesis much more attractive, with the implication 
that the present concept of the conservation of lepton number comes into question. 

Double beta decay [5] in which two neutrinos are emitted (2z/2/3) has been observed 
in a number of cases with lifetimes in the range of 10^^~^^ years. If the Majorana 
hypothesis is correct, it becomes feasible for two virtual neutrinos to annihilate and 
the resulting Oz/2/3 decay would proceed with the two electrons bearing the total decay 
energy [3]. The latter process is expected to proceed at a slower rates with hfetimes 
greater than 10^^ years, with the value depending on the nuclear matrix element and 
the electron neutrino rest mass, as well as the available energy. The observation of this 
process may well be the only way to resolve issues on the basic nature of the neutrino 
and to provide information on the absolute value of its rest mass. 

To extract the mass, it will be crucial to have confidence in the reliability with 
which the nuclear matrix element can be calculated. Most nuclear processes have 
been characterized by theoretical descriptions somewhat empirically, by looking at 
systematics and then developing a model that can describe the systematics consistently. 
This provides a semi-empirical framework for further calculations. For the neutrinoless 
double beta decay, we are unlikely to have this luxury since, even if it is observed, the 
number of cases are likely to be insufficient to establish systematics. It is therefore 
desirable to obtain as much data as possible to characterize the relevant properties of 
the nuclei that are the most likely candidates for the observation of such decays. 

The framework used for estimating the matrix elements has been largely that 
provided by the proton-neutron quasiparticle random phase approximation (QRPA) [5] 
which is successful in describing overall features of nuclear properties. QRPA depends 
on some rather drastic simplifying assumptions and is thus often not quantitatively 
accurate in describing the details of nuclear structure. On the other hand, large model 
spaces can easily be accommodated. Shell- model calculations are also used [B]. The shell 
model can generally describe most features of nuclear structure more quantitatively, 
but there are limits on the range of orbitals that can be included in the calculations. 
This restricted model space requires renormalisation of transition operators and the 
closure approximation [7j is invoked for the virtual intermediate states in neutrinoless 
double beta decay. Detailed aspects of nuclear structure can also be described using so- 
called algebraic approaches, where the interacting boson approximation (IBA) has been 
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found to be very useful. It is therefore not surprising that double beta decay matrix 
elements have also been calculated using the IBA model [8], again utilizing the closure 
approximation. 

From an experimentalists perspective, the dominant question is to understand 
which readily measurable nuclear properties probe aspects of the nuclear Hamiltonian 
that are important for the nuclear matrix elements. These properties can then be 
measured and the results used to investigate the extent to which the theoretical models 
used for the matrix element calculations reproduce them. 

In the past few years, we have carried out a series of experiments designed to 
test some of the nuclear properties that are likely to be relevant to the neutrinoless 
decay process. In this article we review these experiments and the general methods. 
We also attempt to summarize the issues that remain in determining the relevance of 
these properties to the decay process. This article is written from the perspective of 
experimental physicists trying to identify the relevant information. 

2. What Aspects Of Nuclear Structure Matter? 

The question of how easily observable features of nuclei may be relevant to Oz/2/3 decay 
is, somewhat surprisingly, not very well explored. We give some qualitative arguments 
here - but their validity needs to be investigated more rigorously. 

2.1. Relationship to other weak processes. 

The understanding of simple beta decay in nuclei is reasonably complete, especially for 
allowed transitions. The phenomenal constancy of super-allowed transitions between 
mirror nuclei [9] is one of the cornerstones of the field. But decay rates to more 
complicated states, especially in nuclei close to regions of changing deformation, are 
more difficult to predict. 

Matrix elements of the 2z/2/3 decay process are calculated by allowing for the virtual 
excitation of states in the intermediate nucleus between the parent and daughter. Since 
the momentum provided by one of the beta decays is of the order of ~1 MeV/c, only 
the lowest multipolarities (l"*" states) at the lowest excitation energies can contribute 
significantly. If the energies and matrix elements of these intermediate states can be 
measured, for example by charge exchange reactions, a reasonable estimate can be made. 

For neutrinoless Oz/2/3 decay, however, there is a fundamental qualitative and 
quantitative difference. Here the two neutrons decaying into protons emit virtual 
neutrinos that must annihilate in a short distance within the nucleus [TU] and such 
a distance implies high virtual momenta, up to several 100 MeV/c. Consequently, the 
multipolarities involved in the intermediate virtual excitation [TT] can be quite high 
- perhaps up to 8 - and excitation energies in the intermediate nucleus up to ~100 
MeV can play a role in the virtual process. The sensitivity to the structure of the 
intermediate nucleus, which is dominant for the 2i/2/3 process, is not very relevant here 
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[12] . Because of this difference, it does not seem reasonable to assume that a theory 
successfully describing the matrix elements for single or double beta decays with real 
neutrinos represents any assurance that the same method can also reliably predict the 
neutrinoless mode. 

2.2. The relevance of the intermediate virtual excitations. 

Charge-exchange reactions such as (^He,t) can excite particle-hole modes in the 
"giant resonance" region at excitation energies from 10 to 25 MeV. The experimental 
measurements of 1"*" excitations studied in this way have been crucial to calculations of 
the matrix elements for the 2i/2/3 mode that is dominated by the virtual excitation 
of a few 1+ states at low excitation in the intermediate nucleus. For the case of 
the Oz/2/3 process however, this is not the case because of the high virtual momenta 
available in the intermediate vertex. This leads to many higher multipoles and excitation 
energies up to 50-100 MeV, encompassing all relevant particle-hole excitations. For the 
neutrinoless mode, it is not at all clear that results from charge-exchange reactions are 
likely to be useful in constraining the calculations, since all possible modes are within 
easy reach. The question of whether a closure approximation is appropriate has been 
investigated, and it has been found to be satisfied to better than 10% [T^ at least in 
QRPA approaches. 

To the extent that the closure approximation is valid, we can say that the details of 
the intermediate nucleus are not of major relevance in determining the matrix element 
for neutrinoless decay since the response function of nuclei to particle-hole excitations 
in the region relevant to this process has essentially all the possible modes accessible. 

2.3. Initial and final states. 

What nuclear properties can make a difference for the Oi/2/3 decay? The specific 
configurations of the valence nucleons that make up the ground states certainly have to 
make a difference. 

All the candidates for such decays occur between even nuclei with spin-parity 0+. 
The decay involves the destruction of a pair of neutrons and the creation of a proton pair. 
But if there are additional changes to the configuration of nucleons, such as differences 
in the occupancy of the valence orbits greater than just one pair, the re-arrangement of 
many nucleons could inhibit the decay, just as it does for other nuclear processes such as 
single (3 decay or 7 decay, and illustrated schematically in Figure 1. Such effects are not 
within the present framework of QRPA calculations, but they do occur in real nuclei. 
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Figure 1: Caricature of how the rearrangement of nucleons would inhibit transitions. The initial 
state is illustrated on the left, with just two valence orbits. The two upper boxes on the right illustrate 
"allowed" transitions where a pair of neutrons decay into a protons pair in either the same or a different 
orbit. The two lower boxes on the right schematically show configurations such as might arise from a 
rearrangement of valence orbits, changing deformations, or other effects. 

One example would be a change in deformation between the initial state and the 
final one. This involves the rearrangement of many nucleons and is likely to inhibit 
the matrix element for the transition severely, perhaps by orders of magnitude, as it 
does in electromagnetic transitions in the shape isomers in heavy nuclei [13] . Inhibition 
based on the same mechanism can also occur in /3 decay. For example, in the region 
of changing deformation around neutron number = 90, there are large order-of- 
magnitude variations in /t-values for transitions to various 0^ states in the transitional 
nucleus ^^^Sm, depending on whether the decay is from the = 91 1+ state of ^^^Pm, 
or from the = 89 0~ state of ^^^Eu, as is illustrated in Figure 2. 
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Figure 2: The ft-values for /3-decay to the first three 0+ states in the transitional -'^^^Sm nucleus are 
drastically different, depending on whether the decay is from the more deformed N=91 ^^^Eu, or from 
the more spherical N=89 ^^^Pm. For ease of display the ft-values for each decay were divided by the 
sum of the ft values to the three states. 



In this paper we discuss two classes of measurements that we have carried out 
to characterize the initial and final states in the nuclei that are candidates for the 
observation of 0i>2P decays: 

• the populations of various single-particle orbits of valence nucleons, with particular 
attention to the changes in occupancies, and 

• the extent to which correlations between the zero-coupled pairs of like nucleons, 
that lead to initial and final states consistent with the Bardeen-Cooper-Schrieffer 
(BCS) approximation, are confined to the ground states and are similar for the 
parent and daughter nuclei. 

The details of how much specific differences affect the matrix elements are not yet 
clear, but some general considerations will apply. 



2.4- Valence Populations. 

The simplest way to consider differences between initial and final states is to look at the 
decomposition of these states in terms of populations of single-particle orbits of valence 
nucleons. The purpose is, first of all, to check that the changes did not imply some 



Constraining the Oz/2/3 matrix elements by nuclear structure observables 



7 



inhibition through re-arrangement effects and then to identify the specific orbits that 
are changing for both protons and neutrons. 

It has not yet been possible to gain a simple understanding of how specific orbits 
contribute to the process since, in the current QRPA formalism, their population can be 
controlled only indirectly, by modifying the energies of single-particle states for example. 
What has been done [12] suggests that the naive expectations are correct; for example, 
that the largest contributions to the matrix elements occur when the protons are created 
in the same orbits in which the neutrons decayed, for example Ogg/2 to 0g9/2- 

It is qualitatively obvious that the process will be inhibited when there are 
rearrangements in the configurations of nucleons other than the simple change implied 
by the decay, but there is considerable uncertainty in the quantitative estimates of such 
inhibitions. 

2.5. Pair Correlations. 

Both modes of double beta decay involve the conversion of two neutrons into two protons, 
so it may perhaps be expected that correlations between nucleons may be relevant to 
the nuclear matrix elements for decay, as they are to pair-adding or pair-removing 
reactions. One important type of correlation is pairing, where the attractive nucleon- 
nucleon forces lead to the formation of J'^ = 0^ nucleon pairs. At the simplest level, the 
displacement of the even-even and odd-odd mass parabolas due to the pairing energy 
is responsible for the very existence of double beta decay candidates. The attractive 
interactions localize the nucleon pairs spatially and distribute nucleons between single- 
particle orbitals close to the Fermi surface, which thus becomes smeared and orbitals 
acquire partial occupancy. 

If the Oz//5/9 matrix element is written as a sum over the angular momentum J'^ 
of products of neutron pair-annihilation and proton pair-creation operators, the part of 
the matrix element arising from nucleon pairs coupled to total spin 0"^ or J'^ 7^ may 
be separated. Recent calculations [HI [T51 [IS] suggest that while the contributions via 
non-zero pairs, which have a long range, are significant, they tend to have opposite sign 
to the contribution from the 0"*" pairs. There are cancellation effects which diminish 
the long-range components leaving a short-range peak pjj, reminiscent of the wave 
functions of a pair of neutrons in ^H. 

In shell-model treatments, pairing correlations are treated exactly, at least within 
the model space used in the calculation. In QRPA methods, pairing correlations between 
like nucleons are treated separately from other effective interactions via a transformation 
to a quasiparticle regime within the BCS approximation [T7]. However, there are well- 
established circumstances in nuclei where the simple BCS approximation fails, as will be 
discussed in more detail below. How this may alter the calculation of matrix elements 
has not yet been explored. 
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3. Single-Nucleon Transfer and Valence Populations 

In order to compare the differences in valence populations between initial and final 
states, the most suitable techniques are provided by single-nucleon transfer reactions. 
In principle, knockout reactions provide an alternative, but these are not likely to be 
practical for the present purpose because of the energy resolution required. In nucleon- 
adding or nucleon-removing reactions, the angular distributions are characteristic of 
the orbital angular momentum transfer and the cross sections are proportional to 
spectroscopic overlaps between states. These overlaps, or spectroscopic factors, were 
formally defined by Macfarlane and French [18] and obey sum rules that relate their 
sums to the occupancies of single-particle orbits. The summed spectroscopic factors 
for nucleon-adding reactions with a given quantum number reflect the vacancies in the 
corresponding orbit, while the sums for nucleon removing reflect the occupancy. 

3.1. Normalizations and Absolute Spectroscopic Factors 

While there is some uncertainty in the absolute normalization of the reaction 
calculations, it has recently been shown that consistent results can be obtained by 
measuring both nucleon-adding and nucleon-removing reactions on the same target 
[I9l l20l |22] . The method consists of requiring a normalization such that for a given 
orbit characterised by total angular momentum j, the sum of the measured occupancy 
and vacancy on the same target add up to the degeneracy of the orbit 2j + 1. These 
measurements were carried out with special attention to minimizing errors in cross 
sections in order to reduce systematic uncertainties. The neutron-transfer reactions 
were chosen to optimize momentum matching, using (a,^He) reactions for higher £ 
values and (d,p) reactions for lower values. The analyses were carried out within the 
Distorted Wave Born Approximation (DWBA) with the distortions taken from global 
optical potentials. This procedure leads to an independent value of the normalization 
for each target nucleus, and the consistency between four independent determinations 
is illustrated in Figure 3 for the stable Ni isotopes. The four values were then averaged, 
so that the same normalization was used for all the neutron-transfer measurements 
on Ni. This same normalization also gives the neutron occupancies of the valence 
orbits, constituting an independent check of the consistency. These results were not 
very sensitive to the choice of potential parameters in the reaction modelling as long as 
a consistent procedure was used; the uncertainties arising from the choice of potentials 
were at a level of a few percent in the summed spectroscopic factors and under 10% for 
the individual spectroscopic factors. This normalization procedure removes a large part 
of the uncertainty in using transfer reactions to obtain occupancies; it relies on reaction 
theories primarily for the rather small changes in comparing reactions with different 
Q- values. 
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Figure 3: Results from single-neutron transfer reactions on four Ni isotopes. The box on the top 
left illustrates the consistency in the degeneracy, once the normalized i=l spectroscopic factors are 
summed for both occupancies and vacancies. The lower left shows the changing neutron occupancies 
in the valence orbits derived from the transfer data, and their consistency with expectations, with 
the dashed line indicating the expected levels of valence-neutron occupancies. On the top right the 
fragmentation of the £=1 strength is shown, as a function of excitation energy, using negative energies 
for hole states from neutron-removal. The bottom right shows a histogram with the binned strength 
fitted to a Lorentzian, to indicate that not much of the total strength is likely to be missed. 

It has been previously noted, particularly for (e,e'p) reactions, that absolute 
spectroscopic factors fall short of the single-particle value, and this has been understood 
in terms of the limitations of the mean-field approximation arising from high-momentum 
short-range correlations between nucleons [211 [25]. The normalizations obtained from 
the transfer data by the above procedure are similar to the values expected on the basis 
of the quenching of low-lying single-particle strength observed in (e,e'p) experiments. 
The value of this renormalization is apparently a uniform property of nuclear matter, 
and appears to be independent of nucleus, at least to the extent that the limited range 
of target nuclei allows us to make such tests. 

The magnitude of the renormalisation adopted for the occupancy measurement was 
the average of the values determined for individual targets in a given vicinity of nuclei. 
The variation between the values was a few %, as shown in Figure 3. The consistency 
of the normalizations, and the fact that they lead to sensible values of the occupancies 
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to a few tenths of a nucleoli, indicates the level of validity in such a procedure for 
determining changes in valence occupancies. The normalizations depend slightly on the 
optical-model parameters used, but the normalized spectroscopic factors are essentially 
the same for different "reasonable" parameters. 

A great deal of information is available from transfer reaction studies performed 
since the 1960s, which established the foundations of our current understanding of 
nuclei and the single-particle skeleton underlying nuclear structure. Unfortunately, cross 
sections were often not published; in many cases only graphs of angular distributions and 
tables of spectroscopic factors are available. The analyses were carried out with a variety 
of parameters, computer codes, and approximations. Thus the precision and consistency 
needed for quantitative comparisons of data on different nuclei is not readily available 
from the published literature, even where the relevant measurements were carried out 
with great care. 

For the purposes of meaningful comparisons of valence occupancies between initial 
and final states in double beta decay, the older experiments could only be used as a guide 
for identifying i values and to indicate which transitions were strong and which were 
weak. An example of the type of information that has been available is shown in Figure 
4, together with a point indicating the accuracy in absolute cross section obtained in 

m- 




Lab. Angle 



Figure 4: Data taken in an earlier measurement [U] for a transition in ^^Ge(d,p)^^Ge to the 
0.159-McV excited state. The blue points with error bars are values obtained from the figure in this 
reference and their stated absolute errors. The green and black lines represent two DWBA calculations 
with different global distorting parameters. The red bar illustrates the fractional accuracy in the 
measurements that were carried out to help constrain the Oi^2/3 matrix element |19j . 
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3.2. Summary of what has been done, what is under way, and what remains. 

In this section we will review what has been done in determining changes in valence 
populations of nuclei that are candidates for the observation of Oz/2/3 decay. The major 
experiments searching for this exotic decay mode that are currently under way utilize 
^^Ge, ^''°Mo, ^'^°Te, ^^^Xe, and ^^°Nd. Other candidate decays under consideration or 
being tested for suitability are "^^Ca, ^^Se, ^^Zr, ^^^Cd, etc. We have completed the 
measurements on determining the valence-nucleon populations in ™Ge, and are well 
underway for ^°°Mo and ^^''Te systems, with measurements associated with ^^^Xe and 
^^°Nd planned for the near future. In this section we discuss the status of the experiments 
undertaken so far. 



3.2.1. The valence orbit occupancies in '^^Ge and '^^Se. 

The study of this system has been completed and was published in Refs. [121 1^ . 
Some of the general considerations are discussed in more detail for this system in order 
to illustrate the issues involved with all the studies. The germanium and selenium 
isotopes had been studied with (d,p) reactions several decades ago, such that £ values 
and approximate spectroscopic factors were known [21]. Many of the levels that are 
populated in the transfer reactions had also been studied by 7-ray spectroscopy and 
thus spins and parities were confirmed. The active neutron orbits are lp3/2, O/5/2, lpi/2, 
and 05'9/2, with 44 neutrons in ^^Ge. 

For a given £ value, the calculated angular distributions of (d,p) reactions predict 
the peak angles to a reasonable accuracy and the cross sections in Refs. [19] were 
measured with good statistics at these angles. The same procedure was followed for the 
(p,d) reaction, where the bombarding energy was chosen to yield proton and deuteron 
energies as close to those in the (d,p) reaction as possible. To obtain spectroscopic 
factors we rely on the DWBA model primarily to give the energy dependencies. Data 
were obtained for four target nuclei to provide some redundancy in the measurements: 
^^'^^Ge and ''^'^^Se, with 32 and 34 protons for Ge and Se and 42 and 44 neutrons for 
the two isotopes chosen. 

A test similar to that shown in Figure 3 was carried out to give confidence in having 
covered all strengths. Requiring the sum of vacancies and occupancies to be 2j+l yielded 
very similar normalizations: 0.53±.02 for ^=l and 0.54±.06 for £=3. In a chronologically 
later experiment on four nearby Ni isotopes that was already mentioned [22] , comparable 
values were obtained: 0.57±.02 and 0.52±.03 respectively. This consistency is further 
evidence indicating that the procedure followed is valid. The results for the valence 
occupancies in the Ge and Se nuclei are shown in Figure 5. 
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Figure 5 Occupancies for neutron and proton orbitals in ^^Ge and ^^Se measured in Ref.s 120) . 
Not all the spins for the £—1 transitions were known unambiguously and therefore only the total £=1 
strength is shown. 

The proton spectroscopic factors were obtained in a similar fashion, using (^He,d) 
and (d,^He) reactions |20]. However, Ge and Se have 32 and 34 protons respectively, 
just four and six protons beyond Z = 28, and well below Z = 50, in the same major 
shell as the neutrons. Thus the proton-removing (d,^He) reaction is critical for the 
experiment. Since the Z=50 closed shell is far away, the centroid of the transition 
strength for proton addition will be at higher excitation energy (over 3 MeV for 9/2+) 
and tend to be fragmented into many small components. In contrast, for the neutron 
reactions with 42 or 44 neutrons, these nuclei are about midway between shells and 
the neutron-transfer centroids are well below 2 MeV excitation. The normalization 
procedure for proton transfer therefore had to be somewhat different from that followed 
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for neutrons. 

The (d,^He) measurements were done at the RCNP in Osaka with the Grand 
Raiden spectrograph |26]. At this facihty the lowest practical energy was 80 MeV 
for the deuterons, which still provides reasonable momentum matching for the orbits of 
interest. For the (d,^He) reaction, the momentum matched well for i about 2.5, and still 
reasonable for 1 and 4, thus suitable for the transitions to p, f and g orbitals of interest. 
Assuming a single normalization for the three i transfers, one can normalize the DWBA 
calculations by requiring that the occupancies be equal to 4.0 for germanium and 6.0 
for selenium, providing a fourfold redundancy in the normalization. This yielded proton 
occupancies of 3.8, 4.0, 6.1 and 6.2 for the four targets, ^^'^^Ge and ^^'^^Se, indicating 
consistency at the level of around 0.2 nucleons. 

A comparison between the differences in orbital occupancies for neutrons and 
protons from our measurements are shown in Figure 6. The QRPA calculations that 
were done before [27] and after [28] publication of the experimental results are also 
shown in the figure. It is clear that the prior calculations did not fit the measurements, 
but adjusting the assumed single-particle energies in the QRPA considerably improved 
the agreement and the calculated decay rate changed by about a factor of two. 



Differences in Valence Populations '"Ge - '"Se 




Figure 6: Changes in the neutron and proton populations of valence orbits from the measurements 
in Ref.s [Tni[10]j along with those from QRPA calculations from before and after the transfer-reaction 
measurements were published. The changes are shown as positive for both protons and neutrons, even 
though the neutron number decreases. 

3.2.2. The valence orbit occupancies in ^^^Te and ^^^Xe. 

This measurement has been started, but the data analysis is not yet final [21]. The 
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neutron occupancies were determined for i^s.isorpg similar techniques to those 

described in the previous section. Since these nuclei are close to the = 82 neutron 
shell, with 4 and 6 vacancies, the critical measurements are the determinations of the 
neutron vacancies by neutron- adding reactions. The orbits that are important are ld3/2, 
1^5/2, Ohii/2, and possibly 0g7/2. These reactions had been previously studied and the 
strengths seemed to be concentrated in the low-lying states. Because some of the £ 
values were high, the (a,^He) reaction was important to achieve the required momentum 
matching. For the inverse (^He,a) reaction, the i=5 strength was found to be in a single 
state, but the i=2 strength was found to be fragmented and apparently spread to high 
excitation energy, so that it could not be determined reliably. Using a normalization 
deduced only for i=5 and the requirement that the total neutron vacancies add up to 
4 or 6, consistent results were obtained. No evidence for vacancies in the g7/2 orbit was 
seen. 

For the measurement on xenon isotopes, in order to get the thin targets required to 
provide adequate resolution and to avoid the windows associated with a gas-cell target, 
a frozen layer of xenon was formed by spraying gas onto a thin diamond foil held at 
a low temperature [30] • Although this limited the measurement to low beam currents, 
reasonable data were obtained, albeit with somewhat lower resolution than in the Te 
measurement, which used conventional targetry. The summed strengths in both cases 
were consistent with the same normalization. The results on the difference in neutron 
occupancies are shown in Figure 6. 

Experiments to determine the proton occupancies and changes in proton 
occupancies are in the planning stage. This is particularly interesting in view of the 
apparent pair vibration that was observed in (^He,n) studies that will be discussed in 
the Section 4. 



2.5 



Difference in Neutron Occupancy 
I30xe - i30Te 
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Figure 7: Changes in tlie measured neutron populations of valence orbits between ^^°Te and -'^'^"Xe. 
The £—2 strength is probably mostly in the orbit, but not all the spins are known, so this is 
labelled Id. 

3.2.3. The valence orbit occupancies in ^^^Mo and ^^^Ru. 

High resolution measurements of proton-adding and neutron-removing reactions have 
been carried out using the Q3D spectrometer at the Maier-Leibnitz Laboratory in 
Munich |3T]. The reactions were on targets corresponding to the double beta decay 
parent and daughter nuclei, and also on ^°^Ru, which served as a consistency check and 
to help estimate uncertainties better. The (p,d) and (^He,d) reactions were used to 
obtain good momentum matching for low i transfers, and (^He, a) and (a,t) reactions 
used to extract information for higher i values. 

To normalize the spectroscopic factors, the total vacancies for valence protons 
in Mo and Ru were required to be equal to 8.0 and 6.0 respectively, and the total 
occupancy for valence neutrons beyond N=50 to 8.0 for ^°°Mo and ^°^Ru, and 6.0 for 
^°°Ru. The analysis is almost complete and provisional results for the occupancies have 
been presented in Ref. [3T] . 

4. Pairing Correlations and Two-Nucleon Transfer Studies 
4.1. Pairing, Two-Nucleon Transfer, and BCS. 

Since the most important pair correlations in nuclei are those manifest in the short-range 
spatial correlation of 0"*" pairs of identical nucleons, a single-step reaction transferring 
two such nucleons is a specific probe of pairing. Indeed, light-ion induced two-nucleon 
transfer reactions have been used for many years to experimentally access pairing effects. 
Most of the existing data is for the two-neutron-removing (p,t) reaction at bombarding 
energies of several tens of MeV. Considerable information was also obtained for the 
neutron-pair-adding (t,p) reaction in the period from ~1970-90, but this has diminished 
in recent years due to the increased regulations governing the use of radioactive triton 
beams. The analogous light-ion proton-pair-adding reaction (^He,n) is comparatively 
less well studied; the constraints of time-of-flight neutron detection limit the energy 
resolution, and while conditions were optimized in several dedicated facilities, such 
setups are not commonly found today. Not all systems had been studied using (^He,n) 
as carefully as they were for (p,t) and (t,p) reactions. 

The usefulness of two-nucleon transfer in nuclear spectroscopy was realized by 
Yoshida |32], who first analyzed spectroscopic amplitudes for two-neutron transfer 
reactions in a Born approximation. In both or ^He, there is a pair of identical 
s-wave nucleons coupled to 0^. The cross section for pair transfer is enhanced when 
it is to, or from, a state in which BCS correlations result in a coherent mixture of 
pairs. A 0+ pair of correlated nucleons is well localized and thus has good overlap with 
the pair in the mass-3 nucleus. The resulting collective enhancement in the reaction 
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amplitude represents the participation of pairs in the various valence orbitals, adding 
with a common phase. Further, the binding energy of mass-3 nuclei is such that the 
momentum change needed to transfer an i=0 pair is favorable. Simple quantitative 
estimates of the enhanced ground-state yield, compared to a two-quasiparticle excited 
state, yield a^s^gs/ c"gs-s>2qp ~ [321 133] • For a medium mass nucleus, one may estimate 
that the cross section between ground states described by two fully-paired BCS wave 
functions may be enhanced by a factor of around 20-30 compared to the population of 
other non-collective 0^ excited states. 

The gross features of a two-nucleon transfer spectrum should therefore be the strong 
population of the ground state of the final system, with excited 0"*" states populated less 
than a few % of the ground state, assuming that the ground states of target and residual 
nucleus are well described by BCS wave functions. 

4-2. Departures from Simple BCS Behavior, Changing Shapes and Pairing Vibrations. 

The appearance of significant population of excited O'*' states would indicate a departure 
from this simple picture and there are two common situations when this can occur: where 
there are gaps in the single-particle structure that are larger than the pairing interaction 
responsible for the correlations, and where there is a change in the shape of the nuclear 
ground-state. 

Gaps in the single-particle structure close to the Fermi surface gives rise to the 
phenomenon called pairing vibrations, the name arising from their discussion in terms 
of the collectivity associated with fields that promote a pair from one correlated state 
to another |31]. Microscopically, the presence of a gap with an energy larger than that 
associated with the pairing interaction effectively isolates the sets of levels above it from 
being included in the correlated mixture that forms the ground state. Depending on the 
details of structure, these higher single-particle states may then produce another excited, 
pair-correlated 0^ state at higher excitation energy, with significant cross section for 
pair transfer. The effect of the gap is therefore to share the observed BCS pair-transfer 
strength between the ground and the excited 0"*" state. 

Near to significant closed shells such as in ^'^'^Pb, the robust energy gap gives 
rise to well-developed pair-vibrational behaviour and the (p,t) and (t,p) reactions 
populate a near-harmonic set of pairing- vibrational states in neighbouring even-even Pb 
isotopes [35] . Sub-shell gaps in spherical or deformed single-particle energy levels larger 
than the pairing interaction, and/or variations in the strength of the pairing interaction 
among orbits, can also be responsible for pairing-vibration effects. The presence of such 
pairing-vibration departures from a simple BCS behavior indicates that the quantitative 
validity of the simple BCS approximation in the description of the ground-state wave 
functions is altered. This may well have important implications for the calculation of 
matrix elements for double beta decay, perhaps at a level comparable to or larger than 
the splitting of strength seen in two-nucleon transfer. 

In regions of nuclei that span a change in ground-state shape, strong population of 
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excited 0+ states can occur where the overlap between the initial target wave function 
with a final excited-state wave function is large compared to that with the final ground- 
state wave function. Such effects arise in transitional regions, where the ground-state 
shape of the residual nucleus is different from that of the target. In such cases, there is 
often an excited state in the final system whose shape is more similar to that of the initial 
nucleus than the ground state, and thus strong pair transfer is facilitated to the excited 
state. The classic example of this phenomenon is at N=90, and actually involves one of 
the candidates for double beta-decay studies. It will be discussed in more detail below; 
the ill-defined and changing shape will complicate calculations of matrix elements. 

For measurements relevant to double beta decay, two-nucleon transfer experiments 
focus on the population of 0"*" states, starting either with target nuclei that correspond 
to a double beta-decay parent or to a daughter. Excited 0+ final states stand out 
unambiguously because their angular distributions are uniquely forward peaked and the 
cross sections are the measure of spectroscopic overlap. If the final states are more than 
a few percent of the ground-state transition, the pair-removal process in double beta 
decay may also be fragmented and not well described by a theory that fails to take the 
changing shape into account exphcitly. 

4-3. Experimental Information on Pair Transfer in Candidate Nuclei. 

The available information on pair transfer for double beta decay candidates is 
summarised here. 

4.3.1. ^'^Ce-^'^Se 

The nuclear structure of even Ge and Se isotopes close to stabihty is rather well studied 
and evidence for excited 0"*" states has been established. In some of the lighter isotopes, 
two-neutron transfer reactions have shown significant strength populating some of these 
states and this has been interpreted in terms of shape coexistence effects. There are no 
such effects in the A=76 system. 

Recent work [37] was aimed at studies in the vicinity of ^^Ge for double beta decay 
using the (p,t) reaction on targets of '^^'''^Ge and ''^'''^Se. The experiment was focused 
on measuring accurate cross sections at forward angles near the peak of the i = 
angular distribution. Typical spectra are shown in Figure 8. For the targets that are 
relevant to the double beta decay, ^^Ge and its daughter ^^Se, there is no population of 
excited 0"*" states beyond a few percent of the yield of the ground-state transitions. 
The experimental cross sections for the two nuclei, measured for the ground-state 
transitions, are remarkably constant within the 5% experimental uncertainty. These 
features suggest that for the neutrons in ™Ge and ^^Se the quantitative nature of the 
pair correlations does not change appreciably between mother and daughter, and that 
the BCS approximation is reasonable. 
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Excitation Energy (keV) 



Figure 8: Excitation energy spectra for the ^®Ge(p,t) and ^^Se(p,t) reactions from experiments 
described in Ref. ^7\. States with angular distributions consistent with 0+ states are labelled. Peaks 
labelled with an asterisk are target contaminants. There are no excited 0^ states in either spectrum 
at a level of more than a few percent of the ground-state transition. 

There have been no studies of proton-pair transfer on these nuclei. However, 
experiments are currently underway at Notre Dame University to study the (^He,n) 
reaction on targets of '^^'^^Ge p'S j . 

4.3.2. i30reJ30Xe. 

In the ^^°Te-^^°Xe double beta decay system, there is some prior information on 
pair-transfer reactions. The most recent neutron-pair transfer work [32] focused on 
the ^^^'^^°Te(p,t) reactions in the forward angle region that is most sensitive to the 
population of O"*" states. The measurement shows no strong population of excited O"*" 
states, they were all <4% of the yield of the ground state. Measurements were also 
made of the ^^^'^'^°Xe(p,t)^^'''^^^Xe reactions with a frozen xenon target and, while the 
analysis is not yet complete, it seems qualitatively consistent with what was seen with 
tellurium targets. 

However, the situation for proton pairing in tellurium nuclei is very different. The 
(^He,n) proton-pair-adding reaction populates at least one excited O"*" state on targets 
of 122-130'jg -^vith significant strength, ~40% of that of the ground-state transition |10] . 
This striking departure from the BCS expectations is an example of proton pairing 
vibrations. The major shell between Z=50 and 82 has an additional gap, apparently 
separating the first 14 nucleons in the gY/2 and ^5/2 orbitals from the /iii/2, Si/2 and 
c?3/2 with 18 nucleons, creating a subshell at Z = 64. This gap is apparently sufficiently 
large, such that the orbitals above Z = 64 cannot participate in the correlated ground 



Constraining the Oi^2/3 matrix elements by nuclear structure observables 



19 



state and thus do not contribute to the coherent BCS state. This should become more 
apparent in the planned proton occupancy measurements mentioned in Section 3.2.2. 
Such fragmentation of proton BCS strength in these systems raises substantial questions 
regarding the calculations of matrix elements for double beta decay with QRPA, where 
a full BCS state is assumed as a starting point. Shell-model calculations have also been 
carried out for the A = 130 system [6], but it is not clear whether these successfully 
reproduce the observed tellurium proton pair vibrations. 

4.3.3. ^^^Mo-^'^^'Ru 

There is considerable evidence that there is a transition from spherical to deformed 
shapes for nuclei in the A ~ 100 region near = 60 with well-developed ground- 
state rotational bands in A > 100 in zirconium isotopes and A > 102 in molybdenum. 
The effects of this are clearly evident in experimental data on (p,t) and (t,p) reactions 
in the transitional region around the ^'^"Mo double beta decay nucleus and a recent 
measurement has focused on a quantitative comparison of the relevant nuclei |12] . This 
measurement, using (p,t) reactions on targets of ^^'^°°Mo and ^°'^'^'^^Ru, found that 
> 95% of the neutron pair transfer strength is contained in the ground-state transition, 
except for the reaction leading to ^'^Mo. In this case it was observed, in common with 
previous studies, that a state at 735 keV was populated with ~20% of the ground-state 
transition strength. The transitional nature of the ^°°Mo double beta decay parent 
nucleus is likely to present challenges for calculations of the matrix elements. 

Some (^He,n) data exist for stable A ~ 100 nuclei [13]. The data on targets of ^"^Mo 
and ^''^Ru indicated no evidence of excited O"*" states, confirming the simple expectation 
for the proton configurations in the associated ground states. 

4.3.4. ^^^Nd-^^^Sm 

Data for Sm(t,p) reactions [H] showed the first evidence for significant population of 
excited 0^ states in a region characterised by a transition between two different ground- 
state shapes using pair transfer reactions. This is evident in subsequent experiments 
both for neutron-pair adding and removing reactions on samarium isotopes with 
similar phenomena seen in neodymium isotopes [l5l H6] . The data can be interpreted in 
terms of a preference in the reaction to populate a final state with a deformation similar 
to that of the initial state. The effect in pair transfer reactions is most pronounced 
between nuclei with 88 and 90 neutrons; this is identical the change in neutron number 
involved in double beta decay. 
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Figure 9: The figure shows the relative cross section to 0+ final states for the (p,t) reaction on Nd 
and and Sm targets with neutron numbers 88 and 90, normalized to the total strength in each case. 
The removal of a pair of neutrons from N=90 nuclei displays serious fragmentation of the transition 
strength associated with the changing deformation in this step. In the case of pair removal from N=88 
nuclei, connecting two spherical states, the transition proceeds mostly from ground state to ground 
state. The data plotted are taken from the references cited in [3T] . 

The transitional nature of both parent and daughter nuclei in the ^^°Nd-Sm double 
beta decay system is likely to present more serious challenges for the calculations of 
the relevant matrix element. Double beta decay involves the disappearance of a pair of 
neutrons from the parent ground state in a similar fashion to the removal of a neutron 
pair from the target ground state in the (p,t) reaction. If the strength of the (p,t) 
reaction is altered by a factor of ~3 by the presence of deformation, the question arises 
whether there are likely to be drastic changes in the double beta decay process as well. 



4.3.5. i36XeJ365a 

The nuclei associated with the ^^^Xe double beta decay nucleus are less well studied, 
largely due to problems associated with the production of targets. A measurement of 
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the ^^^Ba(p,t) reaction leading to ^^^Ba shows no evidence for a large cross section for 
an excited 0"*" state P?]. 

In proton-pair transfer several transitions were observed with significant strengths 
to O"*" excited states on both ^^^Xe and ^^^Ba, similar to those on ^^°Te discussed above, 
and probably associated with the proton-pair vibration and a subshell at Z=64. As for 
the case of ^^°Te, in this case there are also likely to be problems in the calculation of 
matrix elements within the QRPA formalism. 

5. Conclusions 

The discovery that neutrinos have finite rest mass has led to new interest in the search of 
neutrinoless double beta decay. The developments of new large experiments to search for 
neutrinoless double beta decay with high sensitivity may have increased the probability 
that the observation of the process is on the horizon. The reliability of calculations 
of the associated nuclear matrix elements may soon become a rather crucial issue and 
should be tested in as many ways as possible. 

In this article, two classes of experimental measurements are described where certain 
properties of the nuclear wave functions of parent and daughter nuclei are accessible, 
namely the microscopic decomposition of the occupancy of valence nucleon orbitals and 
the extent to which the correlations between nucleons can be described as the pair 
correlations in the BCS approximation. Both of these aspects would appear to have 
direct influence on the calculation of the matrix elements. 

The experimental methods and the present status of measurements were 
summarized, covering many of the candidates for this type of decay. In at least one 
case where the valence occupancies were measured, calculations did not fit them and, 
when they were adjusted, significant changes in the matrix elements were found. There 
are also several cases where the pair-adding and pair-removing strengths measured by 
transfer reactions are significantly fragmented and therefore not consistent with a simple 
BCS model. No adjustments for such effects have been made in the calculations. 

To understand the implications of these results and to guide future measurements 
some further theoretical insights would be highly desirable. Some of the questions, seen 
from the perspective of experimental measurements, follows: 

• A qualitative understanding of the relative importance of the microscopic differences 
in the configuration of the components of the neutron pairs that decay into proton 
pairs, broken down orbit-by-orbit, would be highly desirable. In other words, are 
the contributions to the matrix element relatively larger where the orbitals are the 
same {e.g. gQ/2 to gg/2 vs. gg/2 to P3/2)? Are they larger for components where 
a zero-node pair decays to a zero-node transition {e.g. O/5/2 to OgQ/2 vs. O/5/2 to 
IP3/2)? The answers to such questions are not clear, yet they are most likely buried 
in the details of calculations that have already been carried out. Such insights could 
be important in guiding which measurements of occupancy are likely to be the most 
important. (There are some calculations in which histograms of contributions of 
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matrix elements are shown [12], but the relative importance also depends on the 
changes in occupancies. Such information must be part of the same calculations, 
but a systematic study, separating the changes in occupancy, is not apparent in the 
literature.) 

• To what extent do the BCS correlations matter? We know that reactions that 
transfer a zero-coupled pair of nucleons the cross sections are related rather simply 
to the pair-creation and pair-destruction operators of pairing (BCS) theory, and 
the coherence among the pairs is critical in determining the cross sections. As was 
discussed, the BCS-like correlations between pairs of nucleons (that is assumed in 
QRPA calculations) is modified substantially in several of the candidate decays and 
this is reflected in the fragmentation of cross sections. Is there a related reduction 
in the matrix elements? 

• Is the closure approximation valid? If not, to what extent does it fail, and what 
aspects of the structure of the intermediate nucleus are relevant to this failure? 
(Closure has been investigated and found to be a good approximation [lOj within 
QRPA. A more universal quantitative estimate of the validity of QRPA would be 
desirable.) 

Answers to these questions would be highly desirable because they would shed light 
on which experimentally accessible nuclear properties the calculations need to be able to 
reproduce in order to help constrain the calculations of matrix elements for neutrinoless 
double beta decay. 

This work was carried supported by the U.K. Science and Technology Facilities 
Council and by the U.S. Department of Energy, Office of Nuclear Physics, under 
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